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Abstract 

We  consider  the  problem  of  preparing  specific  encoded  resource  states  for 
the  toric  code  by  local,  time-independent  interactions  with  a  memoryless 
environment.  We  propose  the  construction  of  such  a  dissipative  encoder  which 
converts  product  states  to  topologically  ordered  ones  while  preserving  logical 
information.  The  corresponding  Liouvillian  is  made  up  of  four  local  Lindblad 
operators.  For  a  qubit  lattice  of  size  LxL,  we  show  that  this  process  prepares 
encoded  states  in  time  O(L),  which  is  optimal.  This  scaling  compares  favorably 
with  known  local  unitary  encoders  for  the  toric  code  which  take  time  of  order 
£2  (L2)  and  require  active  time-dependent  control. 


1.  Introduction  and  main  result 

Dissipation,  while  generally  seen  as  detrimental  for  quantum  computers,  can  nevertheless  be 
a  useful  resource  if  suitably  engineered.  Appropriately  chosen  system-bath  couplings  can 
result  in  non-equilibrium  dynamics  where  initial  states  converge  toward  some  dynamical  steady 
state.  This  kind  of  ‘quantum  reservoir  engineering’  has  been  proposed  as  a  viable  approach  to 
the  experimental  preparation  of  interesting  many -body  states  [11,  18].  Remarkable  examples 
include  the  preparation  of  pure  states  with  long-range  order  in  Bose-Einstein  condensates  [6], 
photonic  arrays  [17],  as  well  as  topologically  ordered  states  [2,  18].  More  generally,  Verstraete 
et  cil  [18]  argued  that,  at  least  in  principle,  an  arbitrary  quantum  computation  can  be  realized  by 
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dissipation.  The  corresponding  process  is  similar  to  Feynman’s  clock  construction  and  has  the 
final  state  of  the  computation  as  its  steady  state.  Subsequent  work  [8]  following  this  program 
proposed  dissipative  gadgets  allowing  the  realization  of  different  dissipative  dynamics  during 
subsequent  time  intervals. 

Here  we  examine  the  dissipative  preparation  of  specific  topologically  ordered  states.  While 
not  realizing  a  fully  dissipative  computation,  this  basic  primitive  could  act  as  a  building  block 
in  a  hybrid  scheme  where  initial  states  for  quantum  computation  are  prepared  by  thermalization 
and  subsequent  computations  are  performed  in  the  usual  framework  of  topological  quantum 
computation.  We  ask  whether  dissipative  processes  can  be  used  to  realize  an  encoder ,  i.e. 
a  map  which  turns  states  on  individual  physical  qubits  into  encoded  (many-qubit)  states.  In 
contrast,  previous  work  only  considered  the  dissipative  preparation  of  some  ground  state  without 
guarantees  on  the  logical  information. 

It  is  worth  mentioning  that  various  unitary  encoders  are  known  for  topologically  ordered 
systems.  For  the  toric  code  [9]  on  an  L  x  L  lattice,  Dennis  el  cil  [5]  gave  a  unitary  circuit  with 
two-local  controlled-not  (CNOT)  gates  of  depth  0(L2)  acting  as  an  encoder.  Bravyi  el  cd  [3] 
showed  for  any  evolution  under  a  local  time-dependent  Hamiltonian  acting  as  an  encoder 
requires  time  at  least  Q(L).  In  turn.  Brown  el  al  [4]  present  a  duality  transformation  from  a 
two-dimensional  (2D)  cluster  state  to  a  topologically  ordered  state,  which  can  be  interpreted  as 
a  geometrically  local  quantum  circuit  of  matching  depth.  Dropping  the  requirement  of  locality, 
Aguado  and  Vidal  [1]  gave  an  encoder  with  depth  D(logL)  with  geometrically  non-local  two- 
qubit  gates. 

The  dissipative  encoder  considered  here  may  be  realized  by  designing  suitable 
system-environment  interactions.  This  is  to  be  contrasted  with  schemes  involving  error 
correction,  which  generally  consist  of  syndrome  extraction  by  measurement  and  associated 
correction  operations.  For  the  toric  code,  an  encoding  procedure  of  this  form  was  given  [7]. 
It  involves  active  error  correction  operations  similar  to  the  minimal  matching  technique  used 
in  [5], 

To  define  the  notion  of  an  encoder  in  more  detail,  consider  a  quantum  error-correcting  code 
Q  =  (C2)®k  C  (C2)0"  encoding  k  logical  qubits  into  n  physical  qubits.  Informally,  an  encoder 
is  a  map  taking  any  state  |'F)  e  (C2)®k  into  its  encoded  version  |40  e  Q  c  (C2)®".  (This  notion 
implicitly  assumes  a  choice  of  basis  of  Q .)  Since  we  are  interested  in  a  physical  system  of  n 
qubits,  we  will  require  the  encoder  to  convert  a  ‘simple’  unencoded  initial  state  into  an  encoded 
logical  state.  That  is,  we  ask  that  for  a  fixed  subset  A\, . . . ,  A*  of  qubits  and  a  fixed  product 
state  | cp)am  <8>  ■  •  •  <8>  \<fi)  a„  on  the  remaining  qubits,  the  encoder  maps 


^  meQ  for  all  I*)  e(C2)0*.  (1) 


We  are  interested  in  encoders  realized  by  evolution  under  a  Markovian  master  equation 


Here  the  Liouvillian  has  Lindblad  form 


j 


with  Lindblad  operators  Lj  acting  locally  on  a  constant  number  of  qubits.  We  ask  whether 
the  completely  positive  trace-preserving  map  (CPTPM)  c,c  generated  by  C  is  an  (approximate) 
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Figure  1.  This  figure  indicates  the  relevant  qubits  in  theorem  1,  and  their  appropriate 
initialization  for  encoding:  qubits  Ay,  A  2  are  initialized  in  the  state  to  be  encoded.  Each 
qubit  in  B  U  B'  is  in  the  state  |0),  while  qubits  in  C  DC'  are  in  the  state  |+).  The  state  of 
the  remaining  qubits  V  can  be  an  arbitrary  (mixed)  state.  We  also  illustrate  the  support 
of  possible  realizations  for  logical  operators  (Xy ,  Zy)  and  (X2,  Z 2)  associated  with  the 
first  and  second  logical  qubit,  respectively.  Our  encoding  procedure  requires  choosing 
realizations  for  each  logical  Pauli  generator  such  that  they  overlap  only  at  the  initial 
unencoded  qubits. 


encoder  for  sufficiently  large  times  t,  i.e.  whether  it  can  realize  the  map  (1).  Our  result  is  the 
following. 

Theorem  1.  Let  Q  C  (C2)®2i~  be  the  toric  code  consisting  of  2 L2  qubits  on  the  edges 
of  a  periodic  L  x  L  lattice.  Consider  the  partition  of  the  qubits  into  disjoint  sets 
AC  B C  B'  U C  U Cr  C  T>  shown  in  figure  1.  That  is 

•  A  —  {/\|,  A  2 }  are  two  neighboring  qubits  having  a  common  adjacent  vertex  u*  and 
plaquette  p *; 

•  B  —  {By,  . . . ,  Bl_  1}  and  C  —  {Cy,  . . . ,  CL_  1}  are  located  along  a  vertical,  respectively 
horizontal,  line  passing  through  Ay; 

•  B'  =  {B[, . . . ,  B'l  x\  and  C  =  {C[, . . . ,  C'L_X)  are  located  along  a  horizontal,  respectively 
vertical,  line  passing  through  A2  and 

•  V  are  the  remaining  2 (L  —  l)2  qubits. 


There  is  a  geometrically  loccd  Liouvillian  C  (with  four- qubit  Lindblad  operators)  such  that  the 
following  holds:  for  any  state  pv  on  V  =  an(j  |vj/}  g  (C2)®2,  and  for  any  e  >  0,  we 

have 


e,£(|^}(^U®|+)(+lSi”1)®|0)(0|J:,(i-1)®po)-|^)(vl/| 


^  € 


whenever 


t  (4 ln(2))L  +  21n(16e-2). 

In  this  expression,  we  use  the  trace  norm  ||A||i  =  tr  V  A M  for  Hermitian  operators. 


(2) 
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We  will  give  a  detailed  description  of  the  relevant  Liouvillian  in  section  2.  The  evolution  dc 
implements  a  continuous-time  version  of  a  local  error  correction  process  somewhat  analogous 
to  Toom’s  rule:  excitations  move  toward  a  single  plaquette/ vertex  where  they  annihilate.  An 
analogous  ground-state  preparation  scheme  for  more  general  frustration-free  Hamiltonians  was 
discussed  in  [18].  However,  guarantees  about  encoded  information  appear  to  be  harder  to  obtain 
in  their  generic  scheme.  Furthermore  because  their  construction  requires  the  injectivity  [14] 
property  of  the  associated  projected  entangled  pair  state  (PEPS)  description  (and  hence  blocking 
of  sites),  the  resulting  locality  of  the  Lindblad  operators  will  be  slightly  worse.  In  contrast, 
our  scheme  directly  exploits  the  stabilizer  structure  of  the  underlying  code,  resulting  in  a 
comparatively  simple  Liouvillian.  Indeed,  our  construction  is  optimal  in  terms  of  locality,  i.e. 
the  number  of  particles  involved  in  each  Lindblad  term  [12].  Related  locality  constraints  for 
pure  steady  states  were  derived4  in  [15,  16].  In  our  setup,  the  initial  product  state  determines  in 
a  transparent  way  which  code  state  is  prepared.  Our  work  goes  beyond  the  mere  characterization 
of  steady  states  by  adding  two  key  ingredients:  the  consideration  of  logical  observables  (which 
are  preserved)  and  an  analysis  of  the  convergence  toward  the  ground  space  of  the  toric  code. 

The  bound  O(L)  on  the  convergence  time  established  by  theorem  1  improves  on  the  0(L2) 
upper  bound  predicted  for  the  analogous  construction  in  [18],  without  a  guarantee  on  the  logical 
information.  In  fact,  it  is  tight:  there  are  initial  states  p  which  thermalize  slowly,  i.e.  c,c{p)  is 
far  away  from  the  code  space  Q  for  any  time  t  <^C  L.  In  a  separate  work  [10],  we  provide  a 
general  no-go  theorem  in  this  direction:  dissipative  state  preparation  of  topologically  ordered 
states  requires  at  least  a  linear  amount  of  time  in  L  if  the  Liouvillian  is  local.  Combined  with 
theorem  1,  this  implies  that  the  construction  presented  here  is  optimal  in  terms  of  preparation 
time  among  the  entire  class  of  local  Liouvillians. 

In  summary,  our  work  shows  that  dissipative  processes  can  be  used  to  implement  an 
encoder  for  the  toric  code.  Intriguingly,  this  encoder  is  more  time-efficient  than  the  best  known 
unitary  circuit.  We  stress,  however,  that  both  types  of  encoders  need  to  be  supplemented  with 
additional  mechanisms  in  the  presence  of  noise,  especially  if  the  encoded  information  is  further 
processed.  As  discussed  in  [13],  local  Liouvillians  such  as  the  one  considered  here  are  not 
suitable  for  the  preservation  of  encoded  information. 

2.  Description  of  the  Liouvillian 

2. 1.  A  Liouvillian  for  a  general  stabilizer  code 

We  first  describe  a  generic  Liouvillian  associated  with  a  stabilizer  code  Q  with  stabilizer 
generators  {Sj}jeS-  We  will  subsequently  specialize  this  to  the  toric  code.  Let  P4  =  \(l  ±  Sj) 
be  the  projections  onto  the  ±1  eigenspaces  of  the  stabilizer  Sj.  The  code  space  Q  is  the  ground 
space  of  the  Hamiltonian 


(the  global  energy  shift  is  introduced  for  convenience).  Our  goal  is  to  implement  a  local  error- 
correction  strategy  by  dissipative  evolution.  Concretely,  we  associate  a  unitary  Pauli  correction 


4  In  terms  of  [16,  definition  2],  our  construction  is  capable  of  preparing  any  state  in  ground  space  of  the  toric  code 
as  a  conditionally  asymptotically  stable  state  of  a  Lindblad  dynamics. 
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operator  C7  with  each  stabilizer  generator  S7  such  that  C7  and  S7  anticommute,  i.e.  [C7 ,  Sj }  =  0. 
We  define  the  CPTPMs 

Tj(p)  =  P*pP*  +  CJPJpPJc]. 

Note  that  by  definition,  7}  lowers  the  energy  of  the  term  Pj  in  the  Hamiltonian,  i.e. 

tr(T/(p)Pj)  ^  tr(pPj)  for  any  state  p.  (3) 

While  after  application  of  Tj,  the  stabilizer  constraint  defined  by  S7  is  satisfied,  its  application 
may  create  a  non-trivial  syndrome  (excitation)  for  a  neighboring  stabilizer  S, &,  k  /  j .  By  design 
(i.e.  the  choice  of  correction  operators  for  the  toric  code  discussed  below),  repeated  application 
of  all  {Tj}jes  eventually  removes  all  excitations,  resulting  in  a  state  supported  on  the  code  space 
Q.  It  will  be  convenient  to  introduce  the  averaged  CPTPM 

TM  =  E  Ti  <4> 

which  randomly  chooses  a  syndrome  and  applies  the  associated  correction  map. 

By  the  correspondence  discussed  in  [20],  each  local  CPTPM  7)  defines  a  local  Liouvillian 
Cj  =  Tj  —  id.  We  are  interested  in  the  evolution  under  C  =  Cj,  which,  by  definition,  is  a 
sum  of  constant-strength  local  terms.  Observe  that  C  =  |<S|(7IV  —  id). 

2.2.  Construction  for  the  toric  code 

We  now  consider  the  toric  code  with  qubits  on  the  edges  of  an  L  x  L  (periodic)  square  lattice. 
We  separate  the  Hamiltonian  into 

H  =  Hip)  +  H(v)  where  H(p>  =  M  L2/  -  J]  Sp  )  and  H(v)  =  ^l  L2I-^  Sv 

\  peS^  /  \  veSiv> 

where  the  former  includes  all  plaquette  and  the  latter  includes  all  vertex  terms.  Here,  we  have 
taken  Sp  and  Sv  to  correspond  to  plaquette  and  vertex  stabilizers,  respectively, 


In  other  words,  Sp  is  the  product  of  Z-type  Pauli  operators  acting  on  the  edges  bounding 
plaquette  p  e  S<p>  and  Sv  is  the  product  of  Pauli  operators  acting  on  the  qubits  incident  to 
vertex  v  e  S(v>.  To  define  the  associated  Pauli  correction  operators  {C,:},:e5(.)  and  {Cp}peS<.Pi,  let 
us  partition  the  set  of  vertices  and  plaquettes  into 

s(p)  =  {p^CS^CStf  and  5(l,)  =  {u*}  U5{B). 

Here  p*  is  a  single  plaquette,  consists  of  the  L  —  1  plaquettes  lying  on  a  fundamental  cycle 
of  the  torus  (which  we  refer  to  as  the  ‘equator’,  running  ‘horizontally’  or  ‘east-west’  along  the 
torus)  on  which  p,t  is  located,  whereas  S(p)  are  the  remaining  L2  —  L  plaquettes.  The  vertex  u* 

as  well  as  the  sets  and  S'11  are  defined  similarly  on  the  dual  lattice,  see  figure  2. 

The  subscript  associated  with  these  sets  indicates  the  direction  of  movement  of  an 
excitation  under  application  of  the  local  correction  map.  For  example,  a  magnetic  excitation 
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(a)  Correction  operations  for  plaquettes 


*  V 

V *  , 

*  > 

^ 

™  v 

(b)  Correction  operations  for  vertices 


Figure  2.  This  figure  defines  the  Liouvillian  C  (see  text).  In  particular,  (a)  defines  the  sets 
SZ?  and  S',1'1 .  The  associated  correction  operation  Cp  consists  of  a  single-qubit  Paufi-X. 
It  moves  magnetic  (plaquette-type)  excitations  to  the  neighboring  plaquette  according 
to  the  indicated  arrow.  Similarly,  (b)  defines  the  sets  and  Electric  (vertex- 
type)  excitations  are  moved  from  one  vertex  to  the  next  according  to  these  arrows.  The 
associated  correction  operator  Cv  is  a  single-qubit  Pauli-Z.  The  qubits  A\Ao  carrying 
the  logical  information  are  indicated  in  red.  They  are  both  part  of  the  special  plaquette 
/;*,  and  incident  to  the  vertex  u*.  No  correction  operation  acts  on  the  qubits  A\A2 . 


(caused  by  an  X -error)  on  a  plaquette  p  G  S',1'1  will  move  to  the  neighboring  plaquette  to  the 


north  of  p  under  application  of  the  correction  map  C; 

,.  That  is,  we  define 

CP  =  \T\ 

for  p  G  S^p> , 

px i 

for  p  G 

c«  =  -fc 

for  v  G  , 

cv= 

for  v  G  S(Z}. 

We  will  set 


r  —  r  —  j 

corresponding  to  a  trivial  correction  operation  (this  is  simply  done  for  convenience).  The 
Liouvillian  £  is  then  defined  as  in  section  2.1. 


3.  Proof  of  theorem  1 


The  proof  of  theorem  1  relies  on  two  basic  statements.  The  first  one  concerns  the  logical 
information  encoded  in  the  state:  for  suitable  initial  states,  this  information  is  preserved  during 
the  evolution.  Recall  that  the  toric  code  has  two  encoded  qubits.  Let 


X\  —  (  Xi>  I  ®  ZAl,  Z\  —  I  (^)  Zc  J  ®  ZA|  ,  Y\—\X\Z\ 

^beB  )  \  ceC 


xA„  y2  =  ix2z2 


*2=  |Q9^I®XA2,  Z2  =  I  QxQ Zc 

\beB'  )  \ceC' 

be  the  two-qubit  logical  Pauli  operators  defined  by  figure  1,  and  let  us  set 
v2  =  {P{P2  |  P,  G  {/,  Xu  Yu  za,  P2  G  {/,  X2,  ?2,  z2}} , 


where  I  is  the  identity  operator.  These  operators  play  a  crucial  role  in  the  following  statement. 
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Lemma  1  (Preservation  of  logical  information).  For  any  two-qubit  Pauli  operator  P  — 
P\  0  Pi,  Pj  e  {/,  Xj,  Yj,  Zj],  let  P  =  P \P2,  P j  e  V2  be  its  corresponding  logical  counterpart. 
Consider  an  initial  state  of  the  form 

Po  =  I'WU®  l+)(+S't_1)  ®  \0)(0\cca~]>  0  Pv,  (5) 

where  the  state  pv  is  arbitrary  on  (C2)®2(X"1)2.  Then 

tr(Pef£(p0))  =  |  P  |  V)  for  all  t  and  all  P  e  V2.  (6) 


Proof.  Let  P  =  P1  <g>  P2  be  arbitrary.  Observe  first  that  the  rhs  of  (6)  is  equal  to 

(V\Pm=tr(Pp0)  (7) 

because  of  the  definition  of  P  (cf  figure  2)  and  the  form  of  p0,  i.e.  the  fact  that  the  qubits  in  1313' 
and  CC  are  +1 -eigenstates  of  the  single-qubit  X  and  Z  operators,  respectively.  In  other  words, 
it  suffices  to  show  the  expectation  value  tr (Pef£(p0))  is  time  independent  for  initial  states  p0  of 
the  form  (5).  We  claim  that  an  even  stronger  statement  holds:  we  have 

(e,£)t (P)  =  P  for  any  P  e  V2,  (8) 


i.e.  any  observables  of  the  form  P  do  not  evolve  in  the  Heisenberg  picture. 

To  prove  (8),  observe  that  the  one -parameter  family  of  unital  maps  {(er£)t}?^ o  is  generated 
by  the  adjoint  £  '  of  the  Liouvillian,  hence  (8)  is  equivalent  to 

C\P)  =  0. 

Since  C  =  J  Cj  is  a  sum  of  Liouvillians  Cj  =  7}  —  id  associated  with  stabilizers  j,  it  suffices 
to  verify  that  for  every  j ,  we  have  T-  ( P )  =  P  or 


p;pp)+p/c;pc/p: 


p. 


(9) 


Because  P  is  a  logical  operator,  it  commutes  with  each  stabilizer,  and  hence  also  the  projections 
P±.  Equality  (9)  is  in  fact  implied  by  [Cj,  P]  =  0,  which  we  can  verify  by  a  case-by-case 
analysis  for  single-qubit  (logical)  operators  P  e  {Xi ,  Z\ ,  X2,  Z2}  (the  general  case  then  follows 
since  a  product  PiP2  commutes  with  Cj  if  each  factor  Pj  does).  Consider  for  example  the  case 
where  P  =  Z\.  We  then  have  to  consider  two  cases: 


(i)  The  correction  j  =  p  is  associated  with  a  plaquette.  In  this  case  Cp  is  a  Pauli-X  (or 
the  identity).  However,  inspection  of  the  support  of  Z\  (see  figure  1)  and  the  location 
of  the  correction  operators  (see  figure  2(a))  reveals  that  no  correction  operation  Cp  acts 
on  the  support  of  Z\,  hence  [Cj,  Z\\  —  0  as  desired. 

(ii)  The  correction  j  =  v  is  associated  with  a  vertex.  In  this  case  Cv  is  a  Pauli-Z  (or  the  identity), 

hence  [Cj,  Zj]  =  0  holds  trivially.  □ 


The  second  fundamental  statement  is  about  convergence  to  the  ground  space. 

Lemma  2  (Convergence  time).  Let  Q  be  the  projection  onto  the  code  space  of  the  toric  code, 
and  let  Q  =  I  —  Q  be  the  projection  onto  the  orthogonal  complement.  Let  p  be  an  arbitrary 
initial  state.  Then  we  have 

tr(Q  V£(p))^e  for  all  t  ^  (41n(2))L  +21n(l/c).  (10) 
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We  have  not  optimized  the  constants  in  this  bound  as  we  are  interested  in  the  overall  (linear) 
scaling  in  L.  The  proof  strategy  is  different  from  the  arguments  in  [18]  and  may  be  of 
independent  interest. 


Proof.  Consider  the  function  /:  <S\{v*,  p*}  —>  N  U  {0}  defined  as  follows: 

for  a  plaquette  p:f(p)  +  1  =  length  of  a  path  moving  north-east  from  p  to  p* 

for  a  vertex  v:  f(v)  +  1  =  length  of  a  path  moving  south-west  from  r  to  r*. 

In  other  words,  the  function  expresses  the  axial  distance  to  v*  and  p*,  respectively.  For  a  vertex 
v  (plaquette  p),  the  quantity  f (v)  ( f(p ))  is  the  number  of  vertices  (plaquettes)  traversed  by 
an  electric  (magnetic)  excitation  on  the  primary  (dual)  lattice  before  reaching  u*  (p*)  along  a 
path  (v  =  Vi,  v2, . . . ,  vf(v),  u*)  (or  (p  =  pu  p2, . . . ,  Pf(P),  p*)).  For  the  special  vertex  v*  and 
the  plaquette  p*,  we  set  f(v*)  =  f  (p*)  =  —  1  for  convenience  (alternatively,  we  could  omit 
the  discussion  of  the  corresponding  trivial  correction  operations  altogether  as  the  stabilizer 
generators  are  linearly  dependent). 

The  key  property  of  the  function  /  is  the  fact  that  it  is  compatible  with  the  way  excitations 
are  propagated  under  the  correction  operations.  More  precisely,  for  each  stabilizer  Sj,  let 
Pred(  j)  be  the  set  of  correction  operations  that  anticommute  with  it,  i.e. 

Pred(j)  :=  {k  G  S  \  k  #  j  and  {Ck,  Sj}+  =  0}. 

Then  /  has  the  property  that 

k  G  Pred(j')  implies  f(k)  ^  /O')  +  1.  (11) 


Namely,  whenever  an  excitation  is  created  by  Ck  at  Sj,  one  can  be  certain  that  a  higher  valued 
excitation  has  been  removed  at  Sk .  We  will  argue  that  for  any  a  /  1 ,  we  have 

tr(QV£(p))  /  e“(1“a"lm)*  ^a/y),  (12) 

jes 


where  m  :  =  max^s  \  {j\k  G  P red (/ ) }  |  is  the  maximal  number  of  stabilizers  a  single  correction 


operator  can  excite.  In  the  case  of  the  toric  code,  m  =  1,  hence  we  obtain  tr(Q  c'£(p))  / 
e~t/222L  (implying  the  claim)  by  choosing  a  =  2  and  observing  that 

L- 1  ,  j  ,  \  2 

=  2a-1  ^  ar+s  =  2a~l 

r,s= 0 


J2a 

jeS 


fU) 


a 


a 


2  a 


2L-1 


< 


(a  -  l)2 


=  2 


2  L 


a=2 


To  prove  (12),  consider  the  observable 

d=  J2  «/0)pJ- 

jeS\{v*,p*} 


Clearly  D  0  GO  for  any  a  ^  1,  hence  it  suffices  to  show  that  the  expectation  value  tr  ( Dc,£  (p)) 
is  upper  bounded  by  the  rhs  of  (12),  i.e. 

tr (DetC(p))  ^  e-(1-“_lw)r  J2afU)-  (13) 

jeS 

Consider  the  Heisenberg  evolution  of  the  projection  operators  {P7 }; .  Since  Pauli  operators  either 
commute  or  anticommute,  a  straightforward  calculation  gives 

Ptrp-,_j0  if[C*,S,]  =  0, 

k  ' ((/-  PJ)P/  -  P  P  if  {Ck,  Sj}  =  0. 
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In  particular,  the  expectation  values  behave  classically  under  the  designed  Liouvillian,  i.e. 
(writing  ( X ),  =  tr (XetC(p))  for  brevity) 


d(P  ~)t  ^ 

^  ((/ -PT)P-),-  <P-Py>( 

ke  Pred(  j) 

<MPJ>*  +  J2  Pk)f  (14) 

kePred(j) 


According  to  (14),  we  have 

^  <-  <£>>, +y>™>  <pn. 

j  kePred(j) 

^  _(D)f+a-i^a/W(P-)f  h 

k  j:kePred(j) 

where  we  used  property  (11)  on  the  function  /.  According  to  the  definition  of  m  and  D,  this 
implies 

d  (D)t  , 

— ^7- —  ^  —(1  —  a~lm)(D)t, 
at 

i.e.  the  expectation  value  decays  exponentially.  The  claim  (13)  then  follows  from  the  fact  that 
P 7  ^  /  for  aH  ./  e  <S  since  these  are  projections,  hence  tr(Dp)  =  ( D)0  ^ 

4  a 


With  lemmas  1  and  2,  we  are  ready  to  prove  our  main  result. 

Proof  of  theorem  1.  Consider  an  initial  state  p0  of  the  form  (5)  and  assume  that  t  ^  (4  ln(2))L  + 
21n(c"1).  By  lemma  2,  we  have  tr(Q"Le'£(po))  ^  £•  With  the  gentle  measurement  lemma  (see 
e.g.  [19,  lemma  9.4.1]),  this  implies 

l|ef£(Po)-p;ili  (15) 

where  p't  is  defined  as  p'  :=  .  Observe  that  the  state  p[  is  supported  entirely  on  the  code 

space  Q. 

Note  that  for  any  two  states  p,  cr  we  have  ||  p  —  a  ||i  =  max||p||^i  tr (P(p  —  <r))  and  therefore 
|  tr(P(p  —  a))  |  ^  ||p  —  a  ||  i  for  any  normalized  operator  P.  Hence  (15)  implies  that  for  any 
normalized  logical  operator  P,  we  have 

|tr(P(ef£(p0)-p;))|^2Vi.  (16) 

Let  |4>)  e  Q  be  the  target  encoded  state,  i.e.  the  state  satisfying  (cf  (7)  and  (8)) 

(ip|P|vp)  =  (\p|p|\p)  —  tr(Ppo)  =  tr (PetC(p0))  for  any  logical  operator  P . 

Combining  this  with  (16)  shows  that  T) (T  and  p[  have  approximately  the  same  expectation 
values  for  logical  operators.  Since  both  states  are  supported  on  the  ground  space  Q  for  which 
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we  have  access  to  a  full  algebra  of  logical  observables  (linear  combinations  of  logical  Pauli 
observables),  one  may  choose5  a  normalized  logical  operator  P  such  that  P Q  =  sgn(|T/)('T|  — 
p[ )  which  achieves  the  maximization  defining  one-norm  distance 

\\mm -p't\U  =  tr[sgn(|*)<*|  -#)(!*>  (*|-pf')]  =  \tr(P(mm-p'))\^2^.  (17) 

Combining  (17)  with  (15)  and  using  the  triangle  inequality,  we  conclude  that 
l|ef£(po)  -  mmu  <  if  f  ^  (4 ln(2))L  +  21n(e_1). 

The  claim  of  theorem  1  follows  immediately  from  this  statement.  □ 

4.  Conclusions  and  outlook 

We  have  presented  a  quasi-local  time-independent  Liouvillian  which  generates  a  dissipative 
encoder:  it  encodes  two  physical  qubits  into  the  ground  space  of  the  toric  code.  Its  key 
features  are  translation-invariance  (up  to  boundary  conditions),  optimal  locality  of  the  Lindblad 
operators,  as  well  as  optimal,  i.e.  linear  convergence  (encoding)  time  of  the  resulting  dissipative 
encoder  as  a  function  of  the  lattice  size.  This  illustrates  the  power  of  engineered  dissipation  for 
state  preparation  in  a  non-trivial  setting  with  degeneracy. 

Since  the  Liouvillian  is  constituted  of  quasi-local  CPTPMs  associated  with  each  stabilizer, 
one  may  interpret6  the  resulting  evolution  clC  as  the  repeated  application  of  the  ‘average’ 
correction  map  (cf  (4)).  This  viewpoint  immediately  results  in  a  preparation  algorithm  for 
surface  code  states.  The  latter  inherits  the  properties  of  the  dissipative  evolution:  it  is  translation- 
invariant  up  to  boundary  conditions,  composed  of  feedback- free  quasi-local  maps  and  requires 
only  a  minimal  number  of  iterations  to  converge.  These  features  make  it  an  attractive  candidate 
for  potential  experimental  realizations. 

4. 1.  Towards  fault  tolerance 

The  most  pressing  open  question  is  to  incorporate  noise  in  the  form  of  control  and  initialization 
imperfections  into  the  design  and  analysis  of  such  an  encoder.  The  holy  grail  is  to  engineer 
a  reasonable,  possibly  time-dependent  Liouvillian  master  equation  such  that  after  a  specified 
initial  period  physical  qubits  are  encoded  with  only  a  constant,  lattice-size  independent 
probability  of  suffering  uncorrectable  errors.  Expecting  further  reduction  seems  impossible 
without  changing  the  assumptions  due  to  the  possibility  of  an  error  occuring  on  the  unencoded 
qubits.  As  it  stands,  our  protocol  L-fold  amplifies  such  an  initial  error  probability  as  certain 
errors  on  the  qubits  BUC  UB'  UC'  are  propagated  into  logical  errors.  Indeed,  we  suspect  that 
a  fully  fault-tolerant  encoding  scheme  would  require  progressively  growing  the  code’s  physical 
support  while  simultaneously  introducing  disssipative  terms  responsible  for  self-correction 
within  this  support.  As  in  [13]  such  terms  would  be  responsible  for  ensuring  that  the  rate  of 
occurence  for  uncorrectable  errors  decreases  with  the  code  size. 

5  If  P  is  the  normalized  local  observable  which  optimally  distinguishes  between  the  unencoded  states  U 1 1  4/ )  ( |  U 
and  U^PfU,  then  P  —  U  ’  P  U  is  the  required  logical  operator. 

6  To  do  so,  simply  expand  the  exponential. 
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